On the relation between entanglement and subsystem Hamiltonians 
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We show that a proportionahty between the entanglement Hamihonian and the Hamiho- 
nian of a subsystem exists near the limit of maximal entanglement under certain conditions. 
Away from that limit, solvable models show that the coupling range differs in both quantities 
and allow to investigate the effect. 



The reduced density matrix (RDM) which describes a subsystem q of a total system in a 
pure quantum state has been the topic of numerous studies. It can be written in the form = 
exp{—'Ha)/Z with an operator Tia which has become known as entanglement Hamiltonian For 
free fermionic or bosonic systems in their ground state, Tia has again free-particle form and can 
be determined explicitly, see [2] for a review. Because of the thermal form of pa, the question 
whether Tia is related to the subsystem Hamiltonian arises naturally. From the exact results, 
one sees that the answer is in general no. For example, segments in non-critical quantum chains 
like the transverse Ising model or a dimerized hopping model lead to a single-particle spectrum in 
Tia which is linear near zero, whereas has energy bands with a gap. Moreover, the low-lying 
eigenfunctions of Ha are concentrated near the boundaries while they are extended in Ha except 
for zero-energy modes. A certain similarity exists only in the critical case, where both spectra 
are asymptotically linear with level spacing 1/lnL and i/L, respectively, where L is the length 
of the subsystem. This allows to define an effective temperature in the RDM js^. However, the 
eigenfunctions of the two Hamiltonians still differ, and also the forms of Ha and Ha in real space, 
see {2]. 

The situation becomes different if the subsystem is translationally invariant, as is the case for 
sublattices in a chain or for a leg of a ladder. Then the eigenfunctions of Ha and Ha are both 
momentum eigenstates and a closer relation is possible, although not necessary. For a transverse 
Ising chain, for example, the sublattices simply decompose into the individual sites, but Ha has 
non-trivial momentum-dependent excitations in the fermionic representation 14]. In other cases, 
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however, correspondences between the spectra were found, see e.g. P, Q] for quantum Hall systems 
and j3] for a Heisenberg ladder. This feature was explained ina recent paper for coupled conformally 
invariant subsystems with left- and right-moving particles 

In the present note, we want to point out that a relation T-L^ ~ can be obtained very simply 
via perturbation theory for a total system formed from two strongly coupled subsystems. This 
is essentially also the case treated in [S]. We also show, for a solvable fermionic system, how 
away from strong subsystem coupling the operator Tia contains longer-range interactions, as found 
numerically in [9] for Heisenberg and AKLT ladders. 

Consider a quantum system made up of two parts with Hamiltonians Hi and H2 coupled via 
the Hamiltonian H' . This could be a ladder with two legs and rungs described by H'. We assume 
H' large and treat Hi + H2 as a perturbation. Then, if |\I'o > is the (non-degenerate) ground state 
of H' , it changes in first order to 

1*0 >= I'J'o > - 2^ l^fe > _^ (1) 

where > the eigenfunctions of H' and Ej. the eigenvalues. We now assume 

(1) There is only coupling to excited states with the same gap A = Ek — Eq 

(2) Both Ha give the same matrix elements, < ^fc|-ffi|^'o >=< ^fc|//2|*o > 
Then 
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1^1 >= 1^-0 > -—^l^fc >< ^-fclFill-o > (2) 
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which can be written 



1^-1 >= > -^Hil^fo > (3) 

where Hi = Hi— < Hi > with < Hi >=< \I'o|i^i|^o >• The total density matrix then is, to first 
order, 

= >< - — \^Hi >< + 1*0 >< *o| ^^ij (4) 
Since Hi operates only in subsystem 1, the trace over subsystem 2 can be taken and leads to 

p\ = Pi - ^{HiPi + PiHi) (5) 

where pi is the RDM for |\I'o >. If now pi is a multiple of the unit matrix, which means that 
l^'o > is maximally entangled, it can be pulled out in front and one can write, exponentiating the 
differences 

p\ = ^eM-^Hi) (6) 
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where Z = tri(l — AHi/ ^) = tri exp(— 4i?i/A). Thus p\ is correctly normahzed to first order and 
one has the relation 

^1 = (7) 

i.e. a direct proportionality between the two Hamiltonians. The quantity A/4 can be viewed as 
an effective temperature and by assumption one is in the high-temperature limit. Condition (2) 
could be weakened to a proportionality between the matrix elements. This would only change the 
prefactor in ([7]). 

The conditions used in the derivation are not as restrictive as they may seem. They are fulfilled, 
for example, for an antiferromagnetic Heisenberg ladder with Hamiltonian 

H = Hi + H2 + H' = jY, S„S„+i + jY, TnTn+i + J'J2 ^^^^ 

n n n 

where the S and T are spin one-half operators. Then the ground state 1^*0 > is a product of singlets 
at the different rungs. Each singlet is maximally entangled and gives a RDM which is 1/2 times 
the 2x2 unit matrix. Each term in Ha has matrix elements to triplet states at two neighbouring 
rungs which leads to A = 2J' . Thus the coupling in T-Li is given hy K = 4J/2J' = 2J/J'. This 
is exactly the result found numerically in [9i] in the limit J <^ J', see Fig. 5(a) there. In their 
notation, K = 2 cos{9) / s'm{9) and one has to consider 9 w tt/2 where K = 2(tt/2 — 6). 

One can ask if the considerations also hold for an anisotropic Heisenberg model. If the rung 
coupling remains isotropic, this is indeed the case, since the singlet-triplet level scheme for each 
rung does not change. However, if H' is of XXZ form, one has two single levels and one doublet. 
Then there are excitations with two different gaps to the spin singlet, which is the lowest state 
throughout the planar region (| < J'^ = J'y = J'). These appear with different pieces of Hi and 
([7]) is changed to 

^i = -^^M, + -^^M (9) 

where = J' + J'z ^'^^ = 2 J'. Thus while the Heisenberg form remains, the anisotropy of T-Li 
is not the same as that of Hi. An exception is the planar case, Hi^z = 0. Then Hi couples only to 
the doublet and the formula ([7]) with the proper gap holds again. This is interesting, because Hi 
is then solvable in terms of fermions whereas the Hamiltonian of the ladder is not. 

n 

The considerations also apply to a fermionic system as treated in [Sj. Consider two species of 
fermions with opposite dispersion and mutual coupling. The Hamiltonian is 

H = Hi + H2 + H' = Y,lq - ^ 7<? b\bq + {a\bq + b\aq) (10) 

q q q 
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where q denotes the momentum. If 7^ = q, this describes two systems with only right- or left- 
moving particles. If 7^ = cosq, it describes a ladder with opposite nearest-neighbour hopping matrix 
elements in the two legs. The coupling term H' is diagonalized by the operators (og it bq) / \f2 and 
gives the two single-particle levels ±5 for each q, thus A = 2b. Moreover, the levels are analogous 
to spin singlets and therefore maximally entangled. The operators have equal matrix elements 
between them. Therefore d?]) holds for large b and Hx is of the form 

^i = ^egajaq (11) 
g 

with Eq = 47g/A = 2jq/b. This is the result found in 

For this system, however, T-Li can be determined exactly and the e follow from the eigenvalues of 
the correlation matrix in the subsystem 1 



121. But because of the translation invariance, this 



matrix is diagonal in momentum space and the eigenvalues are given by the occupation numbers 
Uq =< OgOq >. Diagonalizing (fTOjl with a canonical transformation Og = Uqaq + Vq(3q,bq = —Vqaq + 
UqPq where Uq + Vq = 1, one obtains 



H = Y, u^qialaq - /3t/3g) , Uq = ^7| + (12) 

q 

This gives the occupation numbers 



V 



U^--) (13) 
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and leads to Eg = ln[(l — nq)/nq], or 



= ii±i (14) 

In [8|] this was obtained in a different way. If one considers the other subsystem, Vq is replaced by 
Uq, which changes the sign of Eq but not the RDM spectrum. 

Expanding (fn|) for large b, one reobtains the result Eq = 2^q/b found above. In the opposite 
case, b <^ jg, however, the variation is logarithmic. Eg = 2ln{2^q/b). The variation of Eq with q for 
7q = —cosq is shown in Fig. 1 for several values of b. One sees that the amplitude increases as b 
becomes smaller. At the same time, the curves deviate from a simple cosine function and become 
more rectangular. This is illustrated for b = 0.1 by the dotted line. Near the points q = =t7r/2, one 
is always in the strong-coupling limit and the slope is ±2/(5. 
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FIG. 1: Dispersion relation for the single-particle excitations in Hi for the model ([TU)) with 7^ = — cosq and 
several values of 5. The dotted line shows a cosine for comparison. 
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FIG. 2: Hopping amplitudes corresponding to the dispersion relations in Fig. 1. Shown are the absolute 
values of tn/ti for odd distances n. 

The deviation of Eg from 7^ means that the hopping in is different from that in Hi. In the 
example, where one has nearest-neighbour hopping in Hi, one finds hopping to more distant sites 
in Til- This can be seen directly by expanding ()14p to higher orders. It is more instructive, though, 
to show the result graphically. This is done in Fig. 2, where the amplitudes t„ for hopping to the 
n-th neighbour are plotted, normalized by ti. Shown are only those for odd distances, since the 
other ones are zero. 
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For large 6, one sees a very rapid decrease of the t„ with n and only ti is relevant. As S 
decreases, the decay slows down and also longer-range hopping becomes becomes important. Due 
to the shape of £q, however, the dominant term is always ti. Very similar results were obtained 
numerically in for Heisenberg ladders. Formally, they are connected with the higher orders 
in the perturbation expansion for |^'o >• The effect is reminiscent of the situation for transfer 
matrices in two-dimensional Ising or Gaussian models, where the exact operators in the exponent 
and those of the Hamiltonian limit differ in dispersion relation and coupling range. 

Summing up, we have shown how a proportionality between Tii and Hi can be obtained for 
strongly coupled and maximally entangled subsystems by treating the subsystem Hamiltonians 
in first-order perturbation theory. Whether they describe a critical or a non-critical system does 
not matter, only their smallness enters. The entanglement is decresed only weakly in this case. 
The fermionic example showed explicitly how the situation changes away from the strong-coupling 
limit. In the free-fermion and free-boson case, one can find a number of simple systems, where 
an expression as (jl4p appears. Examples are hopping ladders with alternating rung couplings or 
the BCS model and the (bosonic) Luttinger model as systems of right- and left-moving particles. 
An exception is a homogeneous hopping ladder. There Hi + H2 commutes with H' and does not 
change the wave function. The matrix elements in condition (2) are then of opposite sign. 



Finally, one should mention that formulae similar to ( 1131) . (1141 ) have appeared before in studies 

nn 

of quenches in quantum chains |13l. Il4i|. In this case, one determines the occupation numbers for 
the modes of the new Hamiltonian in the state before the quench, using the appropriate canonical 
transformation. If, for example, one starts in the ground state of a hopping chain with alternating 
site energies ±6 and switches this dimerization off, as done in |l^, the Ug of the modes in the 
homogeneous chain are exactly (|13p with 7^ = cosq. As above, one can then define a thermal 
density matrix and an effective Hamiltonian, and the difference is only that these quantities refer 
to the full system and not to a part of it. The relation Eg = 2jg/6 in this case was already noted 
infl. 
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